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(N ; ABSTRACT 

I We present near-IR narrow -band continuum-subtracted images in the H2 2.122/im, and Br7 2.166//m 

^ . emission lines for a sample of 15 bipolar planetary nebulae. H2 emission was definitely detected for 

C^ I most of the objects in this sample (13 out of 15). The very high H2 detection rate supports the idea that 

r^ ■ bipolar planetary nebulae have important reservoirs of molecular material and offer suitable physical 

conditions for the excitation of H2. The strength of the H2 emission and the H2/Br7 flux ratio are found 
to correlate with the morphology of the bipolar nebulae observed. Bipolar PNe with broad and bright 
0> ! rings exhibit stronger H2 emission than bipolar PNe with narrow twists. High-quality (sub-arcsec) [N ll] 

^ ' and Hq optical images have been used to compare the distribution of the ionized and molecular material. 

^ . The H2 emission lies just outside the optical [N ll] emission zone. 

I . Subject headings: infrared: interstellar medium: lines and bands - interstellar medium: molecules - 

O I planetary nebulae: general 

ti '• 

j^ ! 1. INTRODUCTION able nowadays make the observations in the H2 2.122 ^m 

> Planetary Nebulae (PNe) form when low- and emission line (hereinafter we will refer to this line as H2, 

•O intermediate-mass « 10 Mq) stars evolve off the asymp- ^^^^^P^ ^^en exphcitly mdicated) a powerful tool for the 

^ totic giant branch (AGB) phase to white dwarfs. In this ^^udy of the distribution of molecular gas m PNe. Since 

d short-hved transition, AGB stars eject most of their enve- *e first detection of this hne in the young PN NGC 7027 

■ " " lope which become thick circumstellar shells. The central (^reffers et al. 1976), over 60 other PNe and proto-PNe 

star will later ionize this material and a new PN arises. ^ave been detected m this line (Webster et al. 1988; Zuck- 

While the visible PN is dominated by ionized gas, molec- ^™^" ^ ^atley 1988; Kastner et al. 1994, 1996, and ref- 

ular material may still survive in high-density clumps or ^'"^"'^^^ therein; Hora & Latter 1996), and the number of 

in the photo-dissociation region (PDR) where the optical detections is increasing steadily. Some doubts about the 

depth to ionizing photons is high (Tielens 1993). The ^^l^^^y of the H2 emission hne as a quantitative probe 

presence of residual molecular gas in PNe, remnant of the "^ ^he total molecular mass have been posed (e g. Dm- 

AGB envelope, provides us with important hints about the ^""^^em 1991) as this line can be excited either by fluo- 

mass-loss history and geometry on and at the end of the '"^^'^^"'^^ through the absorption of UV photons from the 

AGB, and offers some insight into the formation mecha- '^^"^ral star in the photo-dissociation front (Black & van 

• ^ p^T Dishoeck 1987; Sternberg & Dalgarno 1989) or by shocks 

Extensive surveys in the CO (2-1) line in PNe (Hug- (^- S- Burton, Hollenbach, & Tielens 1992). The compar- 

gins & Healy 1989; Huggins et al. 1996, and references ^^^^ between CO and H2 emission in PNe (Huggins & 

therein) have shown that molecular envelopes are widely "ealy 1989; Huggins et al. 1996) shows that the H2 emis- 

observable in PNe and constitute a significant component ^^^^ '' ^^^^^^^ ""^l^^ed to the amount of molecular gas in 

of their masses. The large format near-IR arrays avail- 

' Currently at Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 West Green Street, Urbana, IL 61801 
^Currently at Centro Astronomico Hispano-Aleman, Apdo 511, E-04080, Almeria, Spain 
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PNe. 

The CO and H2 observations of PNe have revealed that 
molecular gas is predominantly detected in bipolar PNe, 
although emission is also found among other morpho- 
logical classes. Bipolar PNe have relatively dense equa- 
torial regions which provide a safe haven for molecules 
due to self- and dust-shielding of UV photons. They are 
also associated to high-mass progenitor stars (Peimbert & 
Torres-Peimbert 1983; Zuckerman & Gatley 1988; Cor- 
radi & Schwarz 1995; Huggins et al. 1996) which pre- 
sumably have experienced higher mass ejections. 

The bipolar PN class is not homogeneous as a whole, as 
two distinct types of detailed morphologies exist (Man- 
chado et al. 1996, hereafter M96; Manchado 1997): (1) 
bipolar PNe with a compact core or narrow waist (class 
B); (2) bipolar PNe with a broad ring structure in the waist 
(class Br). Balick (1987) already noticed this morpholog- 
ical difference (the so-called early and late butterfly types) 
and proposed that B PNe evolve into Br PNe. 

Webster et al. (1988) reported observations of a sample 
of 1 1 bipolar PNe and noted that Br PNe have the high- 
est H2 intensity and H2 to Br7 flux ratio. However, their 
results should be used with caution due to poor spatial 
resolution (~ 5" x 14") and sampling. The latter tends to 
maximize the H2 to Br7 flux ratio. This can be seen in 
some individual objects, as NGC7027, whose flux ratio, 
as given by Webster et al., does not account for the H2 
to Br7 large ratio variations over the nebula (Cox et al. 
1997). 

In order to gain more insigth into the correlation be- 
tween H2 emission and the detailed morphology of bipo- 
lar PNe, this paper presents near-IR narrow-band images 
in the H2 and Br7 lines of a sample of 15 bipolar B and Br 
PNe selected from the narrow-band images displayed by 
M96. The H2 and Br7 fluxes of these PNe have been mea- 
sured and accurate H2 intensities and H2 to Br7 flux ratios 
of each subgroup of bipolar PNe have been obtained. The 
results are compared with these previously presented by 
Webster et al. (1988) and the implications on the forma- 
tion and evolution of bipolar PNe are discussed. More- 
over, the comparison of the H2 and optical images from 
the M96 catalog have been used to analyze the relative 
distribution of ionized and molecular material. Finally, a 
detailed description of the individual morphology of some 
peculiar PNe in the sample is given. 

2. OBSERVATIONS 

2.1. Sample selection criteria 

The nebulae were selected according to the following 
criteria: 

1. All the nebulae are bipolar or quadrupolar PNe, 
showing a 'waist' and (at least) two bipolar lobes. 

2. They have dereddened H/3 fluxes in the range 
3.210-13 < Y{up) < 5.110-11 ergcm-^s-i, and 



[N II] morphologies indicative of the presence of 
dense low-ionization equatorial regions. 

3. They have angular sizes larger than 15" to resolve 
them properly, and smaller than 150" to make them 
fit into the field of view. 

2.2. Data acquisition 

Near-IR images of the selected fields were obtained 
on 24 and 25 June 1996 with the MPI fUr Astronomic 
General-Purpose Infrared Camera (MAGIC) attached to 
the f/8 Cassegrain focus of the Calar Alto (CAHA) 2.2m 
telescope. The detector was a Rockwell 256 x 256 pixel^ 
NICM0S3 array. The high resolution imaging mode was 
used with a pixel scale of 0.64" pixel-i , and a field of view 
of 164". The filter parameters are described in Table 1. 

Different observing strategies were followed for com- 
pact and extended sources. Compact sources were ob- 
served rastering them to different locations of the detector, 
so that the image resulting from the combination of dif- 
ferent frames before centering can be used as a sky frame. 
In the case of extended sources, adjacent blank sky re- 
gions observed between consecutive on-source exposures 
were combined to obtain the corresponding sky frame. 
In both cases, the integration time on the object was al- 
ways 100 sec. The images presented here are composed 
of many shorter exposures which were then individually 
sky-subtracted and flat-fielded before centering and aver- 
aging them to produce the final images. The FWHM of 
point sources in the final images, which includes both the 
seeing at A = 2.2;um and instrumental effects, are in the 
range l"l-l"5. The weather was photometric through 
the whole observing run. The sensitivity at 3cr over the 
background was found to be ~ 3.0 10-^ erg cm-^ s~^ sr-i . 

2.3. Data reduction 

As opposed to optical emission from PNe, strongly 
dominated by narrow nebular emission lines, the situa- 
tion is very different in the near-IR, where the contin- 
uum emission sometimes becomes very important and 
can even dominate the observed emission. Unfortunately, 
there were no H2 and Br7 continuum filters available at 
the telescope. Thus, in order to subtract the continuum 
emission from the narrow-band images, we obtained ad- 
ditional broad-band K' images. 

The subtraction of the continuum emission requires 
some detailed explanation. First, we scaled the broad- and 
narrow-band emission using images of the flux standard 
star SA 107-1006 (CasaU & Hawarden 1992) obtained at 
both filters, assuming that the continuum of this star is 
roughly constant through the ^'-band. Then, for each neb- 
ula we subtracted the emission from the H2 and Br7 to the 
continuum image of the nebula. Finally, after scaling, we 
subtracted these non-contaminated by emission lines con- 
tinuum images to the original H2 and Br7 images in or- 
der to obtain continuum-subtracted images in these lines. 
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Tlie H2 and Br7 narrow-band continuum-subtracted im- 
ages for all the PNe are presented in Figure 1. For com- 
parison, the [N II] optical images from M96 are also pre- 
sented at the same scale. 

The contribution from other emission lines to the flux 
observed in the ^'-band has carefully been considered. 
Besides Br7 and H2 2.122 fim, the most prominent emis- 
sion lines in this band are Br8 H I, He I 2.058, and a few 
H2 emission lines at 1.957, 2.033, 2.223 and 2.248 /im. 
After re-scaling the continuum image, the expected con- 
tribution from Br8 H I and He I 2.058 is estimated to be 
~ 2% the intensity of Br7, whereas the contribution of the 
H2 lines is < 5% the intensity of the H2 2.122 fim line. 

Similarly, we have analyzed the contribution of the He I 
2.1 126;um line to the flux observed in our narrow-band H2 
images. Hora, Latter, & Deutsch (1999) presented reli- 
able measurements of the intensity of the He I 2.1126/im 
and Br7 lines for 33 PNe. The average intensity of the 
He I 2.1126;um line is ~ 5% the intensity of the Br7 
line. Therefore, a significant contribution from the He I 
2.1126;um line to our narrow-band H2 images is only ex- 
pected when Br7 is much more intense than H2. Even in 
these cases, the emission in the H2 and He I lines comes 
from different regions in the PN, and therefore it is still 
possible to distinguish between H2 and He I emission. 

The standard star SA 107-1006 was used for flux cal- 
ibration. The measured H2 and Br7 fluxes are given in 
Table 2. Given the inherent high-noise and short-term 
fluctuations of the near-IR sky background, and the un- 
certainties above described affecting the continuum sub- 
traction, the measured fluxes can be affected by relative 
errors which have been estimated to be < 5% for fluxes 
larger than 1.010"^^ ergcm~^s~\ and 5-20% for smaller 
fluxes. 

3. DETECTION OF MOLECULAR HYDROGEN 

Table 2 lists the usual names^ and approximate galac- 
tic coordinates (the PN G name as described by Acker et 
al. 1992) of the objects observed (columns 1 and 2); the 
morphological classification (column 3); the angular size 
of the whole nebula and its ring (columns 4, and 5) in the 
optical as given by M96; and the dereddened H/3 fluxes 
(column 6) as obtained from Acker et al. (1992). Table 2 
also summarizes the main results of this work: the mea- 
sured fluxes in the H2 and Br7 lines (columns 7 and 8); 
the H2 to Br7 flux ratio (column 9); and the H2 surface 
brightness ratio between the lobes and the ring (column 
10), /lr, defined as 
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where Biobe and Bring is a rough estimate of the peak sur- 
face brightness of the lobes and ring, respectively. 

^The acronym 'He' is used through this work to name Henize's objects, despite the fact that 'Hen' is the acronym used in astronomical dictio- 
naries. The former is used in The Strasbourg-ESO Catalogue of Galactic Planetary Nebulae (Acker et al. 1992) and is widely accepted. 



Emission in the H2 images (Fig. 1) was detected above 
3a level in all but one (He 2-428) of the nebulae in our 
sample. The possible contamination due to the He I 
2.1 126 fim line may become relevant only for the objects 
with the lowest H2 to Br7 ratios: M 1-57, M 1-59, M 2- 
46, and He 2-437. In particular, we cannot exclude that 
the emission detected in M 2-46 is due to He I. The de- 
tection of H2 emission is confirmed in M 1-57, M 1-59, 
and He 2-437 according to morphological criteria; the H2 
emission is detected in the lobes of those nebulae, where 
Br7 emission is not detected. In the case of He 2-437, 
where Br7 emission is detected in the lobes, it should be 
noted that the position angle of the emission in the H2 im- 
age differs from that in the Br7 image, but it follows the 
[N II] emission. Since the H2 fluxes of M 1-59 and M 2- 
46 given in table 2 correspond to the whole nebulae, these 
values should be considered as upper limits; the H2 fluxes 
of M 1-57 and He 2-437 were measured in the lobes, ex- 
cluding the central part. Therefore, the emission detected 
in 13 of the bipolar PNe in our sample is connected with 
H2 emission from these nebulae. 

Figure 2-[top] shows the dereddened Br7 fluxes plotted 
against the H/3 fluxes. Adopting case B recombination, 
Te = lO^K and N^ = 10^ cm-^ the theoretical Br7 to H/3 
flux ratio is expected to be ~ 0.03 (Osterbrock 1989). 
This relation is shown as a dashed straight line in Fig. 2- 
[top]. Although a correlation between Br7 and H/3 fluxes 
is found, most of the objects in our sample fall below the 
predicted Br7 to H/3 fluxes ratio; this is probably due to 
the poor accuracy of the H/3 flux determination, based for 
most of the objects on a rough extrapolation from H/3 flux 
determination on a small portion of the object (Acker et al. 
1991) which is likely underestimating the real H/3 flux. 

The measured H2 fluxes are plotted against the Br7 
fluxes in Fig. 2-[bottom\. It is clear that our sample breaks 
down into two different locations on this diagram. In 
agreement with the values of the H2 to Br7 flux ratios 
listed in Table 2, we distinguish two well defined sub- 
groups in our sample: those PNe whose H2 flux is sev- 
eral times (> 5) stronger than the Br7 flux, and those in 
which the Br7 emission clearly dominates. Globally con- 
sidered, the f (H2)/f (Br7) ratios observed span two or- 
ders of magnitude. As we discuss below, these nebulae 
have different morphological and physical properties. 



3. 1. H2 dominated bipolar PNe 

We can include BV 1, K 3-34, K 3-58, K 4-55, M 1- 
75, M 2-52, M 4-14, and M 4-17 in this group. In all 
these cases the emission observed in the H2 line is more 
than 5 times greater than the Br7 emission. The measured 
H2 flux is in the range 5 10"^^ to 4 10~'^erg cm~^ s"^ , with 
peak intensities 1.2-5.3 x 10~^ergcm~^s"^ sr"^ 



Guerrero et al. 



The morphologies observed in the H2 line (Fig. 1) 
closely resemble those shown in the optical [N ll] lines 
(M96). The H2 emission is mainly concentrated in a 
bright ring, the "waist" of the bipolar nebula. The peak 
emission of the ring has a surface brightness more than 
two times brighter than the projected edge of the lobes, as 
listed in Table 2. The surface brightness contrast between 
the ring and lobes is, however, much less extreme than in 
the optical [N ll] A6583 line. In contrast, the much fainter 
Br7 emission is detected only in the central ring, which 
appears sharper in the H2 line than in the Br7 line. 

All of these PNe share some common properties. First, 
all of them exhibit well defined rings in the light of [N ll]. 
Indeed, they all are Br or "quadrupolar" (Q) with ring PNe 
according to the scheme derived from optical observations 
proposed by M96. Note that BV 1 can also be included 
in this group, although its ring is not so evident, as it is 
seen edge-on (Kaler, Chu, & Jacoby 1988; Josselin et al. 
1999). The rings of these PNe are characterized by a high 
surface brightness contrast between the ring itself and its 
innermost part. The emission in the [N ll] and H2 lines 
arises mainly in both cases from the ring, but the center of 
the ring does not show any emission in these lines or very 
little. These PNe have also relatively low density (100 to 
1000 cm~^) (Acker et al. 1992), and the averaged size of 
their rings is 0.22 pc. In addition, their central stars are all 
faint to be detected in the optical images with the excep- 
tion of the very faint central star of M 4-14. 

3.2. Br7 dominated bipolar PNe 

In contrast with the PNe above cited, the bipolar 
PNe He 2-437, M 1-57, M 1-59, NGC 6881, and 
PC 20 are characterized by a small H2 to Br7 flux ratio 
(F(H2)/F(Br7)~ 0.1-0.5). AU of them show emission in 
the Br7 line strongly concentrated in the central, bright 
core. H2 emission is detected both at the central regions 
and along the bipolar lobes in most cases, although the H2 
emission is always fainter than the Br7 emission. 

He 2-437 and M 1-57 show an unresolved core. They 
also exhibit the most extremely low H2 to Br7 flux ratios 
(~ 0.1-0.2), whereas for M 1-59, NGC 6881, and PC 20, 
the H2 to Br7 flux ratios increase up to 0.2-0.5. These 
latter nebulae, classified as Br PNe by M96, have small 
sized central rings, with a small surface brightness con- 
trast between the innermost part of the ring and its edge in 
the [N II] hne. M 1-59, NGC 6881, and PC 20, therefore, 
show intermediate properties between the Br7 emission 
dominated class B PNe with unresolved cores (He 2-437 
and M 1-57) and the H2 emission dominated class Br PNe 
described above. 

The bipolar PNe in this subgroup exhibit different prop- 
erties than the H2 dominated bipolar PNe. Their central 
regions have a high surface brightness, in agreement with 
their larger densities (10^ to 10^ cm"^) (Acker et al. 1992). 
Their averaged size is only 0.076 pc. Finally, many of the 



central stars of this subgroup of PNe are very bright in the 
optical images. 

3.3. Comments on individual objects 

NGC 6881. 

The [N II] image of this nebula (Guerrero & Manchado 
1998) shows a highly coUimated bipolar structure with a 
size of 29f'5 x 5". Each lobe ends as a bright knot where 
the emission in the [N ll] light is enhanced. A second, 
brighter pair of inner lobes extends over 16". Finally, a 
loop-like structure is located 11" towards the southeast 
extreme of the nebula. 

The differences between the [N ll] and H2 images of 
this nebula are obvious. Unlike observed in the optical, a 
pair of very opened lobes is being collimated by a marked 
4" waist which roughly corresponds to the ring observed 
in the optical. The lobes extend 20" to the north-west and 
13" to the south-east, almost coinciding with the south- 
east loop-like structure. As the lobes separate from the 
central waist, their width increases to a maximum dimen- 
sion of 16". As it is illustrated in Figure 3, the corre- 
sponding width of the optical lobes at this position is less 
than 5". Beyond the limit of the northwest lobe, but at the 
same position angle, extended emission is found as a con- 
tinuation of this lobe. The maximum angular distance at 
which we found emission is ^ 50", almost twice the angu- 
lar dimension of the ionized emission. NGC 6881 is also a 
member of the group of peculiar PNe whose H2 morphol- 
ogy does not trace the ionized emission, as it is the case 
of J 900 (Shupe et al. 1995) and NGC 2440 (Latter et al. 
1995; Latter & Hora 1997). 

Ml-57 

H2 emission is marginally detected ^ 30" from the cen- 
tre of the nebula with the same orientation of the 15" sized 
optical lobes (M96). Deeper images, both in the [N ll] and 
H2 lines, are needed to confirm this result. 

He 2-428. 

This is the only bipolar PNe not detected in the H2 line. 
It is an intriguing situation, since it belongs to the Br class 
of M96 and it is similar in many aspects to other nebu- 
lae in this group, whose H2 flux is a few times stronger 
than the Br7 emission. However, while Br7 emission is 
detected in He 2-428, the nebula is not detected in H2. In 
addition, this is the only case of a nebula in our sample 
in which a bright central star is found. Since the central 
stars of nebulae with bright H2 emission in our sample are 
very faint and difficult to detect, it is urged to investigate 
in detail the nature of He 2-428 and its central star. 

4. MOLECULAR AND IONIZED MATERIAL DISTRIBUTION 

In order to make a fair comparison of the distribution 
of the molecular and ionized material in these PNe, their 
[N II] and Ha-i-[N ll] optical images from M96 were 
rebinned to the same pixel scale of the near-IR images 
(0"64 pixel"^). A gaussian filter was used to degrade 
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tlie optical images and make the FWHM of point sources 
in the field similar for both sets of images (If'l- 1''5). 
The contour plot of H2 emission has been overlaid over 
a gray-scale representation of the [N ll] images in Fig. 1- 
[bottom]. 

Most of the H2 morphological features closely trace the 
corresponding [N ll] optical features. The bright conden- 
sations observed in the rings in H2 are also found in the 
[N II] images. Furthermore, the optical design of the lobes 
in [N II] is also outlined in H2. The H2 emission in the 
bipolar lobes is relatively enhanced compared to the opti- 
cal [N II] lines. The /lr ratios given in Table 2 are in the 
range 0.2-0.4; the corresponding ratios in the [N ll] lines 
are ~ 5 times smaller. 

We extracted spatial profiles from the H2 and [N ll] im- 
ages to quantify the dimensions of the molecular and low- 
ionization material distributions. The spatial line profiles, 
extracted along representative directions of the nebulae, 
are shown in Figure 3. They have been used to measure 
the size of well defined features in both the H2 and [N ll] 
images. In particular, the size of the rings measured from 
peak to peak are listed in Table 3. Given the pixel size 
and the expected FWHM for a point sources, the position 
is accurate within O'.'l-O'.'l. The H2 emission from the 
ring closely embraces the [N ll] optical emission in the 
H2 dominated bipolar PNe. Differences in size goes from 
Of'6 to 3f' 9. Furthermore, emission in the H2 extends fur- 
ther than in the [N ll] line. When lobes are mapped, they 
appear to be more prominent in the H2 line {e. g. K 4-55, 
or M 1-75), partly due to the enhancement of H2 emission 
from the lobes described above. A singular case is the H2 
emission detected in the lobes of NGC 6881 described in 
§3.3. In this case, we measured the lobes width at two 
different slit positions. 

Now, it would be desirable to transform the angular 
differences in size to linear sizes. However, the use of 
statistical distances is uncertain in these objects, since its 
particular morphology makes the spherical symmetry as- 
sumption on which those methods are based unreason- 
able (Pottasch 1984; Maciel 1997). As there are no in- 
dividual distance estimates for most of these objects, we 
have adopted the averaged value among those included in 
Acker et al. (1992), and the more recent estimations given 
by Van de Steene & Ziljstra (1994), and Zhang (1995). 
Using these averaged distances, as given in Table 3, the 
thickness of the region located between the [N ll] emis- 
sion peak and the H2 peak (column 7) result to be in a 
narrow range from 7 x 10^^ to 6.4 x 10^^ cm. Thus, the 
H2 emission occurs in a narrow region just outside the 
ionized zone. 

5. DISCUSSION 

We have searched for molecular hydrogen emission in 
a sample of 15 bipolar PNe. The high detection rate 
(~ 90%) found in our sample is the result of its bias to 



bipolar morphologies, since the frequency of detections 
in a morphologically unbiased sample of PNe (e. g. Kast- 
ner et al. 1996) is rather small, ~ 40%. Moreover, almost 
70% of the combined sample of bipolar PNe surveyed by 
Kastner et al. (1996) and this work exhibit H2 emission. 
These numbers support the idea that bipolar PNe are, at 
least partially, ionization-bounded, with important reser- 
voirs of molecular material, and that they offer suitable 
physical conditions for the excitation of H2. 

The H2 emission closely traces the [N ll] ionized emis- 
sion in most of the PNe in our sample. This behavior has 
been reported, although based on mostly 'qualitative' ar- 
guments (Kastner et al. 1996), data of poor (5"- 10") spa- 
tial resolution (Beckwith, Persson, & Gatley 1978; Web- 
ster et al. 1988), or the quantitative analysis of a few large 
PNe, as it is the case of the eUiptical PN NGC 6720 (Zuck- 
erman & Gatley 1988). The improved spatial resolution 
(~ 1") of our narrow -band [N ll] and H2 images has en- 
abled us to make a detailed analysis of the relative distri- 
bution of the ionized and H2 molecular emission in eight 
bipolar PNe. The H2 emission is found in a thin shell 
(< 0.1 times the optical size), just outside the ionized 
[N II] zone both in the rings and the walls of the lobes. 
As pointed out by Kastner et al. (1996), the H2 emission 
is brightest toward the ring of the nebula. Moreover, our 
data reveal that the H2 emission is enhanced in the lobes 
in comparison to the [N ll] optical emission. 

One of the main results derived from this work is that 
a firm connection can now be established between the de- 
tailed nebular morphology of bipolar PNe and the H2 in- 
tensities and H2/Br7 flux ratios observed. The PNe with 
Br morphology show larger H2 to Br7 flux ratios, as well 
as larger H2 intensities. In contrast, B PNe have the small- 
est H2 to Br7 flux ratios and H2 intensities. This re- 
sult strengthens the earlier suggestion by Webster et al. 
(1988) that strong H2 emission is found almost exclu- 
sively among Br PNe. Adding Webster et al.'s sample to 
ours, we find that ~ 80% of the Br PNe surveyed belong 
to the group with strong H2 intensities. 

Why is H2 intensity enhanced in Br PNe? In order to 
assess and clarify the underlying physics of this correla- 
tion, several hypotheses have been considered, although 
all of them have their own share of difficulties. Moreover, 
they do not exclude each other. 

"B and Br PNe are at different evolutionary stages".— 
This idea, originally proposed by Balick (1987), is sup- 
ported by the different physical properties of B and Br 
PNe described in §3.1 and §3.2. The Br PNe smaller den- 
sities and larger sizes suggest larger dynamical ages. In- 
deed, the high surface brightness of the ring of Br PNe 
requires a long time in order to allow the fast wind from 
the central star to have swept up over the ring the material 
that once was inside it. To this effect, the fact that the cen- 
tral stars of Br PNe are faint is consistent with an evolved 
stage. 
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The H2 to Br7 ratio is expected to increase with time 
for PNe from massive progenitors (Natta & HoUenbach 
1998). According to this model, the predicted ratios in- 
crease from ~ 0.1 for young PNe to ~ 10 for old PNe, 
in reasonable agreement with the observed ratios. How- 
ever, Bobrowsky & Zipoy (1989) and Natta & HoUenbach 
(1998) models predict a steep decrease of the H2 intensity 
with time. For massive progenitors (0.696 Mq < Mcore < 
0.836Mq) with a high density molecular shell (> 10^ 
cm"^), Natta & Hollenbach's model predicts the decline 
of the H2 intensity by a factor > 100 in the early nebular 
evolution. The agreement between the theoretical predic- 
tions and the increase of H2 intensity with time is possible 
if a lower initial density of the molecular envelope (< 10^ 
cm~^) is assumed. In this case, it has been shown (Vicini 
et al. 1999) that the intensity of the H2 2.122 fim line does 
not decline for a PN with core mass 0.7 M© as it evolves, 
but it is enhanced (the reader is referred to this paper for 
further details). If the detailed bipolar morphology is in- 
deed related to the evolution of bipolar PNe, the bulk of 
observational data presented by Webster et al. (1988) and 
this paper indicate that low initial densities of the molec- 
ular envelopes are the rule rather than the exception in 
bipolar PNe. 

"The H2 emission in B and Br PNe have different preva- 
lent excitation mechanisms" .- Note that this does not nec- 
essarily exclude the previous evolutionary hypothesis, as 
the prevalent excitation mechanism may be different as 
the PN evolves. 

Shocks may play a significant role in the excitation of 
molecular hydrogen in bipolar PNe. According to Natta & 
HoUenbach (1998), the shock contribution to the intensity 
of the H2 2.122 ^um line is expected to be 



O.lxMRG/RG/lO-^Meyr- 



(2) 



over most of the PN lifetime, where Mrg and /rg are 
the mass-loss rate and filling factor of the red giant wind. 
Therefore, for high mass-loss rates and/or non-isotropic 
red giant winds, shock excitation would dominate the H2 
excitation. Indeed, the H2 2.122 ^um peak intensities pre- 
dicted by shock models (Draine, Roberge, & Dalgarno 
1983; Burton et al. 1992) for the typical expansion veloc- 
ities (30- 100km s~^) of bipolar PNe (Corradi & Schwarz 
1995) are in agreement with the observed peak intensities 
in our sample (1 .2 - 5.3 x 10~^ erg cm"^ s~^ sr"^). 

However, no evidence for shock excitation is found in 
the archetypical Br PN NGC 2346. The prevalent H2 ex- 
citation mechanism in this PN has been shown to be UV 
excitation (Hora et al. 1999; Vicini et al. 1999). It can be 
argued that even with this evidence, the shock excitation 
cannot easily be disregarded, as the presence of fast winds 
and dissociating shocks are able to produce UV photons 
in the shocked wind which would indirectly excite the H2. 

On the other hand, even if shock excitation is important 
in bipolar PNe, it is not clear why H2 intensity should 



be enhanced in Br PNe and not in B PNe. Both sub- 
samples of nebulae exhibit the same kinematical proper- 
ties: the expansion velocity of the lobes of these bipolar 
PNe ranges from 30 to 100 kms~^ whereas rings expand 
more slowly, at 5-25 km s~^ (Guerrero & Manchado, 
unpublished data). Alternatively, the H2 line intensities 
may strongly depend on the density of the material and on 
the shock structure (Shull & HoUenbach 1978; Burton et 
al. 1992), so that the characteristic geometry of Br PNe 
would favor the enhancement of the H2 emission. 

An additional excitation mechanism that must be taken 
into account in the late evolution of the PN is the presence 
of X-rays from the hot central star. Natta & HoUenbach 
(1998), and Vicini et al. (1999) have highhghted the im- 
portance of the X-ray heating in enhancing the predicted 
H2 intensity, specially for PNe with low initial densities of 
the molecular envelope. However, the uncertainty in the 
far UV and soft X-ray spectrum of PNe forces to use this 
result with special caveat. No X-ray emission has been 
detected in bipolar PNe up to now, most likely because of 
the high extinction and low luminosity of the central stars 
of bipolar PNe (Guerrero, Chu, & Gruendl 1999). 

Obviously, a detailed spatially-resolved spectroscopic 
analysis of the H2 emission line spectrum in a large sam- 
ple of bipolar PNe is needed to fully address these issues. 

"B and Br PNe proceed from different evolutionary 
tracks" .- If Br PNe follow a different evolutionary track 
with enhanced mass-loss rates on the AGB, they would 
have much more massive molecular envelopes than B 
PNe. Subsequently, the H2 emission in Br PNe would 
be stronger than in B PNe. The evolution from more mas- 
sive stars than the progenitors of B PNe or from a binary 
system undergoing a common envelope phase may be the 
cause for the enhanced mass-loss rate in Br PNe. The lack 
of observational evidence to support these particular evo- 
lutionary tracks for B and Br PNe would make it a highly 
speculative conclusion. 

6. CONCLUSIONS 

We have presented continuum-subtracted flux- 
calibrated H2 and Br7 images of a sample of 15 bipo- 
lar PNe. The high H2 detection rate within this sample 
strengthens the idea that bipolar PNe are molecular rich 
objects where H2 finds suitable physical conditions for its 
excitation. 

The H2 intensity and H2 to Br7 flux ratio are found to 
change with the detailed morphology of the bipolar PNe; 
PNe with Br morphology exhibit the highest H2 to Br7 
flux ratios, and H2 absolute fluxes. Although we cannot 
discard the possibility of B and Br PNe proceeding from 
different evolutionary tracks, most observational data in- 
dicate that the observed correlation reflects the fact that B 
and Br PNe are in a different evolutionary stage. In order 
to explain the enhancement of the H2 intensity with time, 
the initial density of the molecular envelope must be low. 
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Tlie different H2 intensities and H2 to Br7 flux ratios ob- 
served in each of these two subgroups of bipolar PNe may 
also be connected to different prevalent excitation mech- 
anisms, which can, at the same time, be a consequence of 
the different evolutionary stage. 

The H2 emission from H2-dominated bipolar PNe is 
mainly located in the ring, although fainter H2 emission is 
also found in the lobes. The comparison between H2 and 
[N II] images of these PNe reveals that the H2 emission is 
distributed in a narrow region just outside the ionized ma- 
terial and that it is more enhanced in the lobes compared 
to the ionized [N ll] emission. 
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Table 1 
Filter description 




Filter name 


Central wavelength 
[/um] 


FWHM 

[/im] 


H2 S(l) 1-0 

Br7 

K' 


2.122 
2.166 
2.100 




0.021 
0.022 
0.340 



Table 2 
Properties of the bipolar planetary nebula sample 



Object 


PNG 


Morph. class^ 


PN size^ 


Ring size^ 


F(H/3)2 


F(H2) 


F(Br7) 


F(H2) 
F(Br7) 


/lr 








[arcsec] 


[arcsec] 


[10-14 


erg cm" 


2s-l] 






Ml-57 


022.1-02.4 


B 


24.3 


10.6 


4100 


21 


239 


0.1 


0.14 


Ml-59 


023.9-02.3 


Br 


24.2 


8.9 


5100 


61 


320 


0.2 


0.22 


M2-46 


024.8-02.7 


Q 


30.0 


7.0 


830 


11 


64 


0.2 




PC 20 


031.7 + 01.7 


Br 


11.8 


6.5 




53 


90 


0.6 


0.75 


M4-14 


043.0-03.0 


Q 


27 


8.5 


720 


192 


18 


10.8 


0.30 


He 2-428 


049.4 + 02.4 


Br 


63 


18 


56 




43 






K3-34 


059.0+04.6 


Br 


20.8 


10.5 


32 


47 






0.37 


He 2-437 


061.3 + 03.6 


B 


45 


4.6 


780 


12 


115 


0.1 




Ml-75 


068.8-00.0 


Br 


11.8 


6.5 


230 


389 


68 


5.7 


0.38 


K3-58 


069.6-03.9 


Br 


23 


12.7 


150 


227 


7.2 


32.0 


0.20 


NGC6881 


074.5 + 02.1 


Q 


29.5 


5.0 


2700 


125 


243 


0.5 


0.19 


M4-17 


079.6 + 05.8 


Br 


24 


20 


230 


324 


30 


10.7 


0.28 


K4-55 


084.2 + 01.1 


Br 


71 


9.0 


120 


245 


17 


14.0 


0.29 


M2-52 


103.7 + 00.4 


Br 


60 


23 


3800 


198 


23 


8.6 


0.18 


BVl 


119.3 + 00.3 


Br 


48 


7 


280 


197 









References.— (1) Manchado et al 1996; (2) Acker et al 1992 
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Table 3 
Comparison of the H2 and [N ii] spatial distributions 



Object 


Position angle 


H2 size 


[N II] size 


Relative size H2/[N ll] 


Distance 


Thickness 




[deg] 


[arcsec] 


[arcsec] 




[kpc] 


[10^^ cm] 


Ml-59 


40 


3.2 


3.2 


1.00 


3.0 


<0.3 


M4-14 


140 


6.4 


5.8 


1.11 


5.2 


2.3 


Ml-75 


54 


9.6 


8.3 


1.15 


3.4 


3.3 




135 


17.9 


16.6 


1.08 




3.3 


K3-58 


8 


8.3 


7.0 


1.18 


6.1 


5.9 


NGC 6881 


45^ 


9.6 


4.5 


2.14 


3.1 


11.8 




45'' 


11.5 


4.0^ 


2.90 




17.4 


M4-17 


27 


20.5 


16.6 


1.23 


2.2 


6.4 


K4-55 


4 


7.7 


7.0 


1.09 


1.4 


0.7 




80 


10.9 


9.0 


1.21 




1.8 


M2-52 


62 


7.0 


7.0 


1.00 


3.3 


<0.3 




120 


11.5 


10.9 


1.06 




1.5 



^Offset of 10 arcsec to the south-east. 

''Offset of 10 arcsec to the north-west. 
TWHM of the line profile. 
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Fig. 1 a- Gray-scale representation of the H2 (top), and Br7 (center) continuum-subtracted images, and of the [N ll] emission line 
image (bottom) from the Manchado et al.'s catalog of BV 1 (left), and M 1-75 (rigth). Contour plot of the 2. 122 /im H2 emission is 
overlaid on the [N ll] image whose intensity levels have been set up to allow a fair comparison with the molecular emission. The 
lowest contour is at 3cr. The insets at the corners of the [N ll] images are displayed at low intensity levels to emphasize the fainter 
regions of the nebula. The image scale is the same for all the images of each nebula. North is top, east is left. 



H2 and Br7 emission in bipolar PNe 



11 



o 

CO 

o 

CO 



o 

Q 

> 

.1— I 



80 



60 



40 



20 



—1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 

M 2-52 Ho 



^ 



20 



40 



60 



80 



80 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

~ Br7 


1 




. 


. 


60 


- 


- 


40 


- 


- 


20 


_ 


■ 1 












•^ - 


r\ 


.,.!,,, 


1 



20 40 60 80 




40 



20 



—1 1 1 1 1 1 1 1 1 1 r- 

M 4-14 Ho 



^ 



20 



40 




60 80 20 

Relative RA (arcsec) 



Fig. lb- Same as Fig. la, but for M 2-52 (left), and M4-14 (right). 
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Fig. Ic- Same as Fig. la, but for K4-55 (left), and M4-17 (right). 
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Fig. \d- Same as Fig. la, but for NGC 6881 {left), and M 1-57 {right). 
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Fig. le- Same as Fig. la, but for K 3-58 (left), and K 3-34 (right). 
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Fig. If- Same as Fig. la, but for M 1-59 (left), and PC 20 (right). 



16 



Guerrero et al. 



o 

CO 

o 
CO 



o 

Q 

> 

• r-i 
-t-> 




40 

Relative RA (arcsec) 



Fig. 1^- Same as Fig. la, but for He 2-437 (left), and M 2-46 (right). 



H2 and Br7 emission in bipolar PNe 



17 



40 



o 

m 
o 



o 

P 

> 
-t-> 



20 



1 1 1 1 1 
" He 2-428 


' ' 1 


■ 


■- 


^m 


. 


w 




■i 




^ 


- 


J 




■ 


- 


. c 


- 


mi- m 

1 


1 



20 



40 



20 



40 




"p 


1 1 1 1 1 


■ 


- 


1 1 1 


1 , , . 1 



20 40 

Relative RA (arcsec) 



Fig. Ih- Same as Fig. la, but for He 2-428. 
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Fig. 2- (top) The Br7 fluxes plotted against the H/3 fluxes as given by Acker at al. (1992). The fluxes have been corrected of 
interstellar extinction using the standard extinction law from Pottasch (1984), and the H/? extinction coefficients computed from the 
comparison of the Ha to H/3 flux ratio given in Acker (1992) and the theoretical Balmer decrement (Osterbrock 1989). The errorbar 
associated to the Br7 and H/3 fluxes have been plotted, as well as the predicted theoretical ratio (dashed line), (bottom) The measured 
H2 fluxes plotted against the Br7 fluxes. When no detection was made, an upper limit of the flux is given and indicated by an open 
dot. The sample appears separated in two regions of this diagram characterized by stronger H2 or Br7 fluxes. 
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Fig. 3- Spatial profiles extracted along selected directions comparing the 2.122 ^m H2 (solid line) and [N II] optical (dashed line) 
emission. The position angle and orientations are indicated on each panel. 



